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Abstract 

We report on a search for the radiative decay of T(1S) to the pseudoscalar mesons n and rf 
in (21.2 ± 0.2) x 10 6 T(1S) decays collected with the CLEO III detector at the Cornell Electron 
Storage Ring (CESR). The 77 meson was reconstructed in the three modes 77 — > 77, 77 — > 7r + 7r _ 7r or 
n — > 7r 7r°7T . The rf meson was reconstructed in the mode rf — ► 7r + 7r — 77 with 77 decaying through 
any of the above three modes, and also 77' — > 7/3°, where — > 7T + 7T~. 

Five out of the seven sub-modes are found to be virtually background-free. In four of them we 
find no signal candidates and in one (T(1S) — > 77/', 77' — > 7T + 7r _ ?7, 7/ — > 7T + 7T~7r ) 

there are two good signal candidates, which is insufficient evidence to claim a signal. 

The other two sub-modes (7/ — > 77 and 77' — > 7/j ) are background limited, and show no excess 
of events in their signal regions. We combine the results from different channels and obtain upper 
limits at the 90% C.L. which are B(T(1S) -► 77?) < 1.0 x 10~ 6 and B(T(1S) 77/') < 1.9 x lO" 6 . 
Our limits are an order of magnitude tighter than the previous ones and below the predictions 
made by some theoretical models. 
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I. INTRODUCTION 

The hadronic decays of heavy quarkonia below the threshold for heavy meson pair pro- 
duction are understood to proceed predominantly via three intermediate gluons. One of 
the gluons can be replaced by a photon with a penalty of order the ratio of coupling con- 
stants, a/a s . Such exclusive radiative decays of the heavy vector mesons J/ip and T have 
been the subject of many experimental and theoretical studies. For the experimenter, the 
final states from radiative decays are relatively easy to identify as they have a high energy 
photon, a low multiplicity of other particles, and low background. Theoretically, the ra- 
diative decays of heavy quarkonia into a single light hadron provide a particularly clean 
environment to study the conversion of gluons into hadrons, and thus their study is a di- 
rect test of QCD. T(1S) — > 77/ is one such candidate channel. This decay channel has 
been observed to be produced in the J /%b charmonium system (the l 3 Si state of cc) with 
B(J/ip — > 777') = (4.71 ± 0.27) x 1CT 3 [1]. Naive scaling predicts that decay rates for ra- 
diative T(1S) decays are suppressed by the factor (q b m c /q c m b ) 2 rs 1/40 with respect to the 
corresponding J/ip radiative decays. This factor arises because the quark-photon coupling 
is proportional to the electric charge, and the quark propagator is roughly proportional 
to 1/m for low momentum quarks. Taking into account the total widths [l| of J /if) and 
T(1S), the branching fraction of a particular T(1S) radiative decay mode is expected to be 
around 0.04 of the corresponding J /if) branching fraction. However, the CLEO search [2| for 
T(1S) — > 77/' in 61.3 pb -1 of data collected with the CLEO II detector found no signal in 
this mode, and resulted in a 90% confidence level upper limit of 1.6 x 10 -5 for the branching 
fraction T(1S) — > 777', an order of magnitude smaller than this expectation. 

The two-body decay T(1S) — >■ 7/2(1270) has been observed |3j in the older CLEO II 
T(1S) analysis, and this observation has been confirmed 0, @], with much greater statistics, 
in CLEO III data. The measurement <B(T(1S) -> 7/ 2 (1270)) = (10.2 ±1.0) x 10" 5 , from the 
combination of the two CLEO III measurements, is 0.074±0.010 times the corresponding J/if> 
decay mode, showing a deviation of roughly a factor of two from the naive scaling estimates. 
In radiative J/ip decays the ratio of rf to ^(1270) production is 3.4 ± 0.4. If the same 
ratio held in T(1S), the 7/ channel would be clearly visible. The channel T(1S) — > 77/ has 
received significant theoretical attention. This channel has been observed in J /if) decays [l| 
with the branching fraction of (9.8 ± 1.0) x 10~ 4 , a value smaller by a factor of five than 
B(J/ip — > 777'). The previous CLEO search of Y(1S) decays produced an upper limit of 
2.1 x 10 -5 at the 90% confidence level for this mode [6|. 

Several authors have tried to explain the lack of signals in radiative T(1S) decays into 
pseudoscalar mesons, using a variety of models which produce branching fraction predic- 
tions of 10~ 6 to 10 -4 . Employing the Vector Meson Dominance Model (VDM), Intemann [3] 
predicts the branching fractions for the heavy vector meson radiative decay into light pseu- 
doscalar mesons. Using the mixing mechanism of 7/, rf with the as-yet-unobserved pseu- 
doscalar resonance rj b , Chao [§[ first calculated the mixing angle X VVb in order to estimate 
the radiative branching fractions. Baier and Grozin [9| showed that for light vector mesons 
(such as J ftp) there might be an additional "anomaly" diagram that contributes significantly 
to the radiative decays. Noting that VDM has no direct relation to QCD as the fundamental 
theory of strong interactions, and referring to [7j, Ma tries to address the problem by using 
factorization at tree level with NRQCD matrix elements to describe the heavy vector meson 
portion multiplied by a set of twist- 2 and twist-3 gluonic distribution amplitudes pTol ] . 
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II. DETECTOR AND DATA SAMPLE 



This study is based upon data collected by the CLEO III detector at the Cornell Electron 
Storage Ring (CESR). CLEO III is a versatile multi-purpose particle detector described fully 
elsewhere . Centered on the e + e~ interaction region of CESR, the inner detector consists 
of a silicon strip vertex detector and a wire drift chamber measuring the momentum vectors 
and the ionization energy losses (dE/dx) of charged tracks based on their trajectories in 
the presence of a 1.5T solenoidal magnetic field. The silicon vertex detector and the drift 
chamber tracking system together achieve a charged particle momentum resolution of 0.35% 
(1%) at 1 GeV/c (5 GeV/c) and a fractional dE/dx resolution of 6% for hadrons and 5% for 
electrons. Beyond the drift chamber is a Ring Imaging Cherenkov Detector, RICH, which 
covers 80% of the solid angle and is used to further identify charged particles by giving for 
each mass hypothesis the fit likelihood to the measured Cherenkov radiation pattern. After 
the RICH is a Csl crystal calorimeter that covers 93% of the solid angle, allowing both 
photon detection and electron suppression. The calorimeter provides an energy resolution 
of 2.2% (1.5%) for 1 GeV (5 GeV) photons. Beyond the calorimeter is a superconducting 
solenoidal coil providing the magnetic field, followed by iron flux return plates with wire 
chambers interspersed at 3, 5, and 7 hadronic interaction lengths (at normal incidence) to 
provide muon identification. 

The data sample has an integrated luminosity of 1.13 fb" 1 taken at the T(1S) en ergy 
y/s = 9.46 GeV, which corresponds to iVr(is) = 21.2 ± 0.2 million T(1S) decays [L2j. 
The efficiencies for decay chain reconstruction were obtained from Monte Carlo simulated 
radiative events generated with the (1 + cos 2 6) angular distribution expected for decays 
T(1S) — > 7 + pseudoscalar. The Monte Carlo simulation of the detector response was based 
upon GEANT [l3j], an d simulation events were processed in an identical fashion to data. 

III. EVENT SELECTION AND RESULTS 

In our search for T(1S) — > 777 and T(1S) — > 77/, we reconstruct 77 mesons in the modes 
V ~~ * 77; V ~^ 7T + 7r~7r°, and 77 — > 7r 7r°7r°; the latter two will collectively be referred to as 
77 — > 37T. We reconstruct the 7/ meson in the mode 777T + 7T~ with 77 decaying in any of the 
above modes, and in addition, the mode 7/ — > 7p°, where p° — > tt + tt~. From the CLEO II 
studies [2], @] we expected five out of the seven modes under investigation to be relatively 
background free and so we employ minimal selection criteria to maximize sensitivity and 
minimize possible systematic biases. The other two, 77 — > 77 and 7/ — > 7/) , have large 
branching fractions, but also large backgrounds, and so our event selection for these modes 
aims to decrease the background with a corresponding loss of efficiency 

Our general analysis strategy is to reconstruct the complete decay chain ensuring that 
none of the constituent tracks or showers have been used more than once, then kinematically 
constrain the intermediate tt° and 77 meson candidates to their nominal masses and finally 
require the event to be consistent with having the 4-momentum of the initial e + e~ system. 
Multiply-reconstructed T(1S) candidates in an event, a problem of varying severity from 
mode to mode, is dealt with by selecting the combination with lowest Xibtai' the sum °f 
chi-squared of the 4-momentum constraint (xfm) an d chi-squared of all the mass-constraints 
involved in a particular decay chain. For example, there are four mass-constraints involved 
in the decay chain T(1S) — > 777'; 77 — > 7t 7t 7t°, three tt° mass-constraints and one rj mass- 
constraint. The mode T(1S) — > 777; 77 — > 7r°7r 7r° is an exception in which we preferred to 
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accept the 77 — > 7r 7r°7r° candidate having the lowest S% = J2i S% i( with S n> i = (m 77 — m n o)/ a 77 
of the zth 7r°candidate. The yield is obtained by counting the number of final state 77 
or 7/ candidates within our acceptance mass window defined as the invariant mass region 
centered around the mean value and providing 98% signal acceptance as determined from 
signal Monte Carlo. Whenever possible, an event vertex is calculated using the information 
from the charged tracks, and the 4-momentum of the photon candidates is then recalculated, 
assuming that the showers originate from the event vertex rather than the origin of the CLEO 
coordinate system. This produces an improvement in the 77 and 7/ candidates' invariant mass 
resolution of roughly 10%, leading to a slight increase in the sensitivity of the measurement. 

The CLEO III trigger [lij relies upon two components: (1) the tracking-based "axial" 
and "stereo" triggers derived from the signals on the 16 axial layers of the drift chamber, 
and the signals registered on the chamber's 31 stereo layers, and (2) the calorimeter-based 
trigger derived from the energy deposition in the Csl crystal calorimeter. The events for 
the "all neutral" modes T(1S) — > 777; 77 — > 77 and T(1S) — > 777; 77 — > 7r°7r°7r are collected by 
the calorimeter-based trigger condition requiring two high energy back-to-back showers. We 
demand that triggered events meet the following analysis requirements: (a) a high energy 
calorimeter shower not associated with a charged track, having a lateral profile consistent 
with being a photon, and having a measured energy greater than 4.0 GeV must be present; 
(b) there must be the correct number of pairs of oppositely charged, good quality tracks 
with usable dE/dx information. The efficiency of these requirements is more than 60% in 
modes involving charged tracks and approximately 54% and 45% for cases where 77 — > 77 
and 77 — > 3n°, respectively. 

The photon candidates we use in forming tt° and rj — > 77 candidates have minimum 
energy depositions of 30 MeV and 50 MeV, respectively. All photon candidates are required 
to be not associated to charged tracks, and at least one of the photon candidates of each 
pair must have a lateral profile consistent with that expected for a photon. The photon 
candidates we use in reconstructing the 77 meson in the 77 mode must be detected either 
in the fiducial barrel or the fiducial endcap 1 calorimeter region only. These candidates are 
then kinematically constrained to the nominal meson mass, the exception being T(1S) — > 
777; T) —>■ 77, where no mass-constraining was done to the 77 candidate, because we examine 
m 77 in this mode to determine our yield. 

The 77 candidates in the mode 7r + 7r~7r° are built by first forcing pairs of oppositely charged 
quality tracks to originate from a common vertex. The ir° candidate having invariant mass 
within 7<t 77 is then added to complete the reconstruction of 77 — > 7r + 7r~7r° candidates. The 
charged tracks are required to be consistent with being pions by adding the pion hypothesis 
SdE/dx = (dE / dx (measured) — dE / dx (expected)) / o dE/dx m quadrature for two tracks and 
requiring the sum of S dE ^ dx to be less than 16. 

In the case of 77 — > 7r 7r°7r°, the 77 candidate is simply built by adding three different 
7T° candidates, where no constituent photon candidate contributes more than once in a 
candidate 77 — > 7r°7r 7r° reconstruction. The 7r° candidates are selected by requiring S n < 10.0. 
In order to increase the efficiency in this mode, an exception was made to the fiducial region 
requirement, and photons in the gap between the barrel and endcap fiducial regions were 



1 The fiducial regions of the barrel and endcap are defined by | cos(0)| < 0.78 and 0.85 < | cos(0)| < 0.95, 
respectively; the region between the barrel fiducial region and the endcap fiducial region is not used due 
to its relatively poor resolution. 
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FIG. 1: Candidate r\ — ► 7r + 7r~7r° (top) and r\ — ► 7r 7r 7r° (bottom) invariant mass distributions from 
T(1S) data. The large number of events near 780MeV/c 2 (top) is due to the abundant process 
e + e~ — ► 7W. No events are observed in our acceptance region, bounded by the arrows. 
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A. The Decay T — > 777, 77 — > 3tt 

The T candidate in the mode 777 is formed by combining a high-energy photon (E > 
4 GeV) with the 77 candidate, requiring that this photon is not a daughter of the 77 candidate. 
The T candidate is then subjected to the 4-momentum constraint of the initial e + e~ system. 
In the case of 77 — >• 37r, multiply reconstructed T candidates were restricted by selecting only 
one candidate. For r\ — > 7r + 7r~7r°, we select the candidate with the lowest XTotai; the sum 
of chi-squared of the 4-momentum constraint and chi-squared of the mass-constraint to 
the 7T° candidate. For 77 — > 7r°7r 7r°, we select the candidate with the smallest S^. The 
selected T candidate is further required to satisfy the 4-momentum consistency criterion, 
restricting Xp4 < 100 for rj — > 7r + 7r~7r° and a less stringent cut of 200 for rj — > 7r 7r 7r° 
measurements. In addition, we limit the number of reconstructed calorimeter showers for 
the mode T(1S) — > 777; 77 — > 7r 7r 7r° to minimize backgrounds such as e + e~ 70 where 
— > K$Kl without jeopardizing the signal efficiency. 

From Monte Carlo simulations, the overall reconstruction efficiencies, e,, for each channel 
are determined to be (28.5 ± 4.3)% and (11.8 ± 1.9)% for the decay chains T — > jr),T) — ► 
7T + 7r _ 7r° and T — > 'yi],i] — ► 7r 7r 7r , respectively. The uncertainties in the efficiency include 
the Monte Carlo samples' statistical uncertainty and our estimate of possible systematic 
biases, which are discussed further in Section HVl 

We find no candidate events within our acceptance invariant mass window for the search 
T(1S) — > 777, i] — ► 37T. The invariant mass distributions for candidate 77 — > 7r + 7r~7r° and 
i] — > 7T 7r 7r°, after imposing all the selection criteria are shown in Figure [TJ 
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B. The Decay T — > 777, 77 — > 77 

The 3-photon final state resulting from T(1S) — > 777; 77 — > 77 is dominated by the QED 
process e + e~ — > 777. Our selection criteria of loosely reconstructing an 77 —>■ 77 meson 
and requiring the x 2 of 4-momentum constraint on the T(1S) meson formed by adding a 
hard-photon to be < 200 are not sufficient to suppress this background. The QED back- 
ground, however has a distinct feature - the two photons having energies E^ and E\ used in 
reconstructing the 77 candidate have a large energy asymmetry, where asymmetry is defined 
as (Ehi — Eio) I (Ehi + Eio) . Real 77 mesons are expected to have an approximately uniform 
distribution of asymmetry in the range (0,1). We require the asymmetry to be less than 
0.8. To further discriminate between the signal and the background, we used a neural net 
approach. 

The input to the neural net is a vector of six variables, namely the measured energy 
and the polar angle 9 of each of the three calorimeter showers used in the reconstruction 
chain. The training sample is comprised of 20,000 simulated signal and background events 
in equal proportion. The simulated e + e~ — > 777 background events have a high-energy 
photon (E > 4 GeV), 77 invariant mass for the two lower-energy photons in the range 
0.4-0.7 GeV/c 2 , and energy asymmetry less than 0.8. 

For our final selection, we choose neural-net output with 51% efficiency while rejecting 
86% of the background. The combined efficiency of our selection criteria for this mode is 
(23.8 ± 2.4)%, which includes possible systematic biases and statistical uncertainties from 
the simulation. The resulting 77 invariant mass distribution from T(1S) data is fit, as shown 
in Figure (2J to a double Gaussian function, whose mass and widths are fixed to values found 
from signal Monte Carlo data, along with a second order polynomial background function. 
From this likelihood fit, we obtain —2.3 ±8.7 events; consistent with zero. We then perform 
the same likelihood fit multiple times fixing the signal area to different values, assigning 
each of the fits a probability proportional to e _x / 2 , where x 2 is obtained from the likelihood 
fit. The resulting probability distribution is normalized and numerically integrated up to 
90% of the area to obtain the yield at 90% confidence level. Our limit thus obtained is 14.5 
events at 90% confidence level. 

C. The Decay T — > 77/, 7/ — > r]iT + 7r~ 

Reconstruction of the decay chains T(1S) — > 777', where 7/ — > ?77r + 7r~, builds on the 
search T(1S) — > 7r/The reconstructed 77 candidate is constrained to the nominal 77 mass. 
The mass-constrained 77 candidate is further combined with a pair of oppositely charged 
quality tracks by forcing the tracks and the 77 candidate to originate from a common vertex. 
In reconstruction of r/; 77 — > 7r + 7r~7r°, care is exercised to ensure that no track is used more 
than once in the decay chain. The high energy photon is combined with the 77' candidate 
to build an T candidate which is further constrained to the 4-momentum of the initial 
e + e~ system. In the reconstruction chain 77'; 77 — > 77, the T candidate with the lowest sum 
of chi-squared to the 4-momentum constraint (xp 4 ) combined with the chi-squared of the 
mass-constraint to the 77 candidate (x 2 ) is accepted as the representative T candidate in the 
reconstructed event. In the modes where 77 — > 3tt, the 7r° mass-constraint chi-squared, x 2 o, 
also contributes to the Xibtai- 

To ensure that only good quality 77 candidates participate in the decay chain, the x% 
values of "77 — > all neutral" candidates are required to be less than 200. Owing to the better 
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FIG. 2: Invariant mass distribution of 77 candidates in T(1S) data for the mode T(1S) — ► 777; 77 — > 
77, overlaid with fits using a) floating area (solid red) yielding — 2.3 ± 8.7 events, and b) area fixed 
to 14.5 events (dashed blue), the upper limit corresponding to 90% C.L. 

measurements of charged track momenta, this criterion is more stringent {x^< 100) in the 
case of 77 — ► 7r + 7r~7r°. The targeted efficiency (around 99%) of this requirement is achieved 
in all three cases. 

The charged tracks used in reconstructing 77' candidates have to be consistent with the 
pion hypothesis. We again require the sum of squared S^E/dx added in quadrature to be less 
than 16 for both the two track and four track cases. The efficiency of this requirement alone 
is around 99%. 

The selected T candidate is further required to satisfy the 4-momentum consistency 
criterion, restricting %p 4 < 100 in the 77 — > 77 case and a less stringent value of 200 for 
77 — ► 37T. The overall reconstruction efficiencies of our selection criteria as determined from 
signal Monte Carlo simulations are (35.3 ± 5.2)%, (24.5 ± 2.2)% and (14.4 ± 2.9)% for 77 
decays to 77, 7r + 7r~7r° and 37r°, respectively. 

After these selection criteria, we find no candidate events in the modes T(1S) — > 777'; 77 — ► 
77 and T(1S) — > jTf',T] — > 7r 7r 7r°, as shown in Figure O However, in the mode T(1S) — > 
jrf] 77 — > 7T + 7r _ 7r , we find two good candidate events passing our selection criteria as shown 
in Figure [3j These two events have been looked at in detail and appear to be good signal 
events. However, they are insufficient to allow us to claim a positive signal, as no candidate 
events are observed in the modes T(1S) — > 777'; 77 — > 77 and T(1S) — > 777'; 77 — > 7r 7r 7r , each 
providing higher sensitivity than the decay chain T(1S) — > 777'; 77 — > 7r + 7r _ 7r°. 



D. The Decay T — > 777', 77' — ► 7/j 

The reconstruction scheme for the decay chain T(1S) — > 777'; 77' — > 7p° is slightly different 
from those previously described. We first build p° candidates by forcing pairs of oppositely 
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FIG. 3: Invariant mass distributions of rjTT + ir^ candidates from T(1S) data. The rj candidate is 
constrained to the nominal rj meson mass. No events are observed in the signal box for rj — ► 77 
(top) and rj — > 7T 7r 7r° (bottom); two signal events are observed for 77 — ► 7r + 7r~7r (middle). 



charged tracks to originate from a common vertex. Next, we add a photon candidate (which 
we refer to as the "soft shower" having energy E s in contrast with the high energy radiative 
photon) not associated with charged tracks, and having a lateral profile consistent with 
being a photon, to build rf candidates. To obtain the maximum yield, we neither restrict 
the energy E s of the photon nor the invariant mass of the p° candidate at this stage. A 
high energy photon is then added, ensuring that the soft shower and high energy photon are 
distinct, to build the T candidate. The T candidate is then constrained to the 4-momentum 
of the initial e + e~ system and the candidate with the lowest Xp4 value is selected. 

The candidate r\' invariant mass resolution is vastly improved due to the mass-constraints 
on the candidate ir° and 77 mesons in 77' — > r]7i + 7i^ decays. In reconstruction of 7/ — > 7/9 , 
a significant improvement in candidate 7/ invariant mass resolution (« 30%) as well as the 
energy resolution of the soft shower is achieved by performing the 4-momentum constraint 
on the T candidate. 

Particle identification in the channel 7/ — > 7p° is achieved by demanding the combined 
RICH and dE/dx likelihood for the pion hypothesis be greater than the combined likelihood 
for each of the electron, kaon and proton hypotheses. Copiously produced QED processes 
such as e + e~ — > r yye + e~ are suppressed by imposing an electron veto, requiring that \E/p — 
1.0 1 > 0.05, where p is the measured momentum and E is the associated calorimeter energy of 
the charged track. QED events of the type e + e~ — > 77/i + /i _ are suppressed by requiring that 
neither track registers a hit five hadronic interaction lengths deep into the muon detector 
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FIG. 4: Invariant mass distribution of 7p° candidates in T(1S) data for the mode T(1S) — ► 777'; 
77' — ► 7/3° overlaid with fits using a) floating area (solid red) yielding —3.1 db 5.3 events, and b) area 
fixed to 8.6 events (dashed blue), corresponding to the upper limit at 90% C.L. 



system. Continuum background of the type e + e~ — > 77P is suppressed by demanding 
E a > 100 MeV. Finally, the event is ensured to be complete by demanding %p 4 < 100. The 
overall efficiency of the selection criteria for this mode is (40.1 ± 2.1)%, including possible 
systematic uncertainties and the statistical uncertainty of the Monte Carlo sample. 

Although highly efficient, our selection criteria are not sufficient to suppress the smooth 
continuum background from the reaction e + e~ — * 77P . The candidate rj' — > 7p° invariant 
mass distribution after our selection criteria, shown in Figure HJ is fit to a double Gaussian 
function over a floating polynomial background function of order one. The parameters of the 
double Gaussian function are fixed to the values obtained from a fit to signal Monte Carlo 
and the area is left to float. The likelihood fit yields —3.1 ± 5.3 events, which is consistent 
with zero. In the absence of a clear signal, we determine the upper limit yield as we do in 
the case of T(1S) — ► 777; 77 — > 77, and find an upper limit at 90% confidence level of 8.6 
events. 



IV. SYSTEMATIC UNCERTAINTIES AND COMBINED UPPER LIMITS 

Since we do not have a signal in any of the modes, and since the kinematic efficiency is 
near-maximal, statistical uncertainties dominate over systematic uncertainties. By compar- 
ison of the expected yield of the QED process e + e~ — > 777 with the calculated cross-section 
for this process, we estimate the uncertainty on the trigger simulation for "all neutral" 
modes to be 4.5%. For modes with only two charged tracks, we have studied the QED pro- 
cesses e + e~ — > 7p° and e + e~ — > 70, and assign a 13% uncertainty on the efficiency due to 
possible trigger mismodeling. For events with many charged tracks, we assign a systematic 
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TABLE I: Contributions to systematic uncertainties in the efficiencies for T(1S) — * jrj' (upper 
half) and T(1S) — > jr/ (lower half). The uncertainties are expressed as relative percentages and 
combined in quadrature. 



Uncertainty source 


j/;jj-> 77 


77' ; ?7 — > 7T+7T TT° 




rf — * 7p° 


Trigger mismodeling 


13 


1 


13 


1 


Track reconstruction 


2 


4 


2 


2 


Calorimeter response 


5 


5 


15 


2.5 


Analysis cuts 


4 


5.7 


4 


3.9 


Monte Carlo statistics 


1.0 


1.6 


2.4 


1.0 


Combined uncertainty 


14.7 


8.8 


20.4 


5.2 


Uncertainty source 


v 77 


T] — ► 7T^7r 7T^ 


T) 7r 7r°7r° 




Trigger mismodeling 


4.5 


13 


4.5 




Track reconstruction 




2 






Calorimeter response 


5 


5 


15 




Analysis cuts 


7 


4 






Monte Carlo statistics 


1.3 


1.2 


1.7 




Combined uncertainty 


9.8 


15.2 


16.0 





uncertainty of 1% as the relevant trigger lines are very well understood, redundant, and 
very efficient. We assign 1% uncertainty per track in charged track reconstruction based 



upon CLEO studies [15j of low-multiplicity events, and 2.5% systematic uncertainty per 
photon from mismodeling of calorimeter response which translates to 5% uncertainty per 
meson (vt° and 77) decaying into 77, again based upon CLEO studies [l5|. The systematic 
uncertainty in SdE/dx f° r two tracks added in quadrature (as in T(1S) — > 777; 77 — > 7T + 7T~7r ) 
was evaluated to be 4% by considering the efficiency difference of this requirement in Monte 
Carlo and data samples of e + e~ — > jou. Consequently, we assign 4% and 5.7% uncertainty to 
the reconstruction efficiencies of modes involving two and four charged tracks, respectively, 
excepting 7/ — > 7/9° where this requirement was not imposed. For the mode 7/ — > 7p°, the 
systematic uncertainty in the efficiency of analysis cuts, found to be 3.9%, was evaluated by 
comparing the efficiency difference in Monte Carlo and data by studying the p° signal due to 
the QED processes. For the neural-net cut in the mode T(1S) — > 777; 77 — > 77, we studied the 
efficiency in QED e + e _ — > 777 simulated events and the real data dominated by the same 
QED process for a wide range of neural-net output values. We find a maximum difference of 
7% in these two numbers, which we take as a conservative estimate of the associated system- 
atic uncertainty. The systematic uncertainties for various rj and 7/ decay modes are listed in 
Table [H These uncertainties were added in quadrature, along with the statistical error due 
to the limited size of Monte Carlo samples, to obtain the overall systematic uncertainties in 
the efficiencies. 

The systematic uncertainties in efficiencies, uncertainties in the product branching ratios, 
and the statistical uncertainty in the number of T(1S) decays, Nrns), are incorporated [l6| 
by a "toy" Monte Carlo procedure to obtain smeared likelihood distributions for the branch- 
ing fraction in each mode, £>(Y(1S) — > 7P) = iVp/(ej • Bpj ■ Nr(is)), where P = 77,77', and 
6j and Bp t i denote the efficiency and branching fractions of the ith mode. To obtain the 
smeared likelihood distribution £p,i, the experiment is performed multiple times, randomly 
selecting A^p from the likelihood function appropriate for each mode 2 and then dividing by 



2 For modes with zero or few observed events, the appropriate likelihood function is generated from Poisson 
statistics. For the background limited modes r\ — > 77 and 77' — > 7/? , we already have the likelihood 



ff 



TABLE II: Results of the search for T(1S) — * jr/' and T(1S) — * 777. Results include statistical and 
systematic uncertainties, as described in the text. The combined limit is obtained after including 
the systematic uncertainties. 







v';v -» 77 


rf\ 1] — > 7T+7T 7T 


V '; V ->Tr TT TT 


n ' - 7P o 


Observed events 

77 j 2 

Reconstruction efficiency (%) 
S(T(1S) -» yrf)(90% C.L.) a 
S(T(1S) -> 7»)')(90% C.L.) b 






17.5 ± 0.6 
35.2 ± 5.2 

< 1.8 x 1CT 6 

< 1.9 x 10" 6 


2 

10.0 ± 0.4 
24.5 ± 2.2 

< 10.3 X 10- 6 

< 10.4 X 10~ 6 




14.4 ± 0.5 
14.4 ± 2.9 

< 5.2 X 10~ 6 

< 5.8 X 10- 6 


-3.1 ±5.3 
29.5 ± 1.0 
40.1 ± 2.1 

< 3.4 x 10~ 6 

< 3.4 x 10~ 6 


Combined limit on B(T(1S) - 


. 7V) 




< 1.9 x 10- 


-b 








T) -» 77 


T] — ► n~^n~ 7T W 


v _> w o^o 




Observed events 
Bn,i% 

Reconstruction efficiency (%) 
B(T(1S) -> 7r;)(90% C.L.)° 
B(T(1S) -» 7r;)(90% C.L.) 6 




-2.3 ±8.7 
39.4 ±0.3 
23.8 ±2.4 

< 7.3 x 10~ 6 

< 7.4 x 10" 6 




22.6 ± 0.4 
28.5 ± 2.9 

< 1.7 X 10~ 6 

< 1.8 X 10" 6 




32.5 ±0.3 
11.8 ± 1.9 

< 2.8 x 10~ 6 

< 2.9 x 10" 6 




Combined limit on B(T(1S) - 


' IV) 




< 1.0 x lO -13 







a excluding systematic uncertainties 
^including systematic uncertainties 



0990307-005 



t — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 








'_ combined 


(b) \ 


/ yv'- ^ 


r\' ,r\^n + K-K° = - 




r j 

n',ri— >7r. 7r°7r.° _ 


r \ \\ 

1 i i i 1 1 \ 1 1 1 1 \ 1 


-../ i 

1 1 1 1 1 1 



20 



10 



15 



20 



Branching Fraction (10~ 6 ) 



FIG. 5: Likelihood distributions as a function of branching fraction for the decay mode T(1S) — > 777 
(left) and T(1S) — > 777' (right). All distributions are smeared by respective systematic uncertain- 
ties and normalized to the same area. The solid black curve denotes the combined likelihood 
distribution. 

the sensitivity factor q • Bp^ ■ iVx(is), where each term is picked from a Gaussian distribution 
about their mean values with the appropriate standard deviation. 



function which we used in calculating the upper limit of the observed number of events at 90% CL. 



12 



The combined likelihood distribution for £>(Y(1S) — > 7P) is derived as Cp = Y\i^v,i 
which is summed up to 90% of the area in the physically allowed region to obtain the upper 
limit branching fraction for T(1S) — > 7P. From the constituent Cp t i and the combined Cp 
as shown in Figure 0, we obtain upper limits on £>(Y(1S) — > 777) of 7.4 x 10~ 6 , 1.8 x 1CT 6 , 
2.9 x 10~ 6 , and 1.0 x 10~ 6 for 77 decaying into 77, 7r + 7r~7r°, 7r°7r 7r°, and all three combined, 
respectively. We obtain upper limits for £>(Y(1S) — > 777') of 1.9 x 10~ 6 , 10.4 x 10~ 6 , 5.8 x 10~ 6 , 
and 3.4 x 10~ 6 for 77 decaying into 77, tt + tt~tt°, 7r 7r 7r°, and 77' — > 7p°, respectively. The 
combined upper limit for /3(Y(1S) — > 777') is 1.9 x 10 -6 , a value larger than one of the 
sub-modes (T(1S) — > 777'; 77 — > 77), due to the two candidate events in T(1S) — > 77/; 77 — > 
7r + 7r _ 7r . The numbers of observed events, detection efficiencies and upper limits are listed 
in Table ED 

V. SUMMARY AND CONCLUSION 

We report on a new search for the radiative decay of T(1S) to the pseudoscalar mesons 77 
and 77' in 21.2 x 10 6 T(1S) decays collected with the CLEO III detector. The 77 meson was 
reconstructed in the three modes 77 — > 77, 77 — > 7r + 7r~7r or 77 — > 7r 7T 7r . The 77' meson was 
reconstructed either in the mode rj' —>■ 7p° or v( — > tt + 7t _ ?7 with 77 decaying through any of 
the above three modes. All these modes except for 7/ — > 7p° had earlier been investigated 
in CLEO II data amounting to iV-r(is) = 1-45 x 10 6 T(1S) mesons and resulted in previous 
upper limits B(T(1S) -> 777') < 1.6 x 10~ 5 and B(T(1S) -> 77/) < 2.1 x 10" 5 at 90% C.L. 
These limits were already smaller than the naive predictions based upon the scaling of the 
decay rate for the corresponding JM radiative decay mode by the factor (qt„m c /q c mb) 2 , and 
also the model of Korner et at, [17], whose perturbative QCD approach predictions for 
B(J/ip — ► jX) where X = 77,77', f% as well as /3(T(1S) — > 7/2) agree with experimental 
results. 

With a CLEO III data sample 14.6 times as large as the CLEO II data sample, we find no 
convincing signal in any of the modes. Based purely upon the luminosities, we would expect 
the new upper limits to be scaled down by a factor of between 14.6 (in background-free 
modes) and \/14.6 in background dominated modes if the two CLEO detectors (CLEO II 
and CLEO III) offered similar particle detection efficiencies. In the search for T(1S) — ► 77/ 
we find no hint of a signal, and manage to reduce the limit by an even larger factor. In 
the search for T(1S) — ► 777', however, we find two clean candidate events in the channel 
T(1S) — ► 777'; 77 — > 7r + 7T~7r°, which, though we cannot claim them as signal, do indicate the 
possibility that we are close to the sensitivity necessary to obtain a positive result. Because 
of these two events, our combined limit for T (IS) — > 777' is not reduced by as large a factor as 
the luminosity ratio, and in fact is looser than that which would be obtained if we analyzed 
one sub-mode (T(1S) — > 777/577 — > 77) alone. In this analysis we found upper limits which 
we report at 90% confidence level as 

B(T(1S) ^777) < 1.0 x 10~ 6 , 

B(T(1S) ^77/) < 1.9 x 10~ 6 . 

Our results are sensitive enough to test the appropriateness of the pseudoscalar mixing 
approach as pursued by Chao [s| , where mixing angles among various pseudoscalars including 
rjb are calculated. Then, using a calculation for the Ml transition T — ► 777^, he predicts 
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B(T(1S) -> 777) = 1 x 1(T 6 and B(T(1S) -> 777') = 6 x 1(T 5 . Our limit for T(1S) -> 7V is 
significantly smaller than Chao's prediction and does not support his approach. 

The sensitivity challenge posed by both the extended vector dominance model and the 
higher twist approach of Ma are beyond our reach. In extended VDM, Intemann predicts 
1.3 x 1(T 7 < B(T(1S) -> 777) < 6.3 x 1(T 7 and 5.3 x 1(T 7 < B(T(1S) -> 77/) < 2-5 x KT 6 , 
where the two limits are determined by having either destructive or constructive interference, 
respectively, between the terms involving T(1S) and T(2S). Even if it is determined that 
the amplitudes are added constructively, our limit remains higher than the VDM prediction 
for T(1S) — > 777. 

Ma's prediction of £>(Y(1S) — > 777') ~ 1.7 x 10~ 6 is consistent with our result. However, 
his prediction for £>(Y(1S) — > 777) « 3.3 x 1(C 7 is a factor of ~ 3 smaller than our limit. 
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